Simple model of 1/f noise 



B. Kaulakys 

Institute of Theoretical Physics and Astronomy, A. Gostauto 12, 2600 

Vilnius, Lithuania 
(Received 1 March 1998) 



Abstract 

Simple analytically solvable model of 1/f noise is proposed. The model con- 
sists of one or few particles moving in the closed contour. The drift period of 
the particle round the contour fluctuates about some average value, e.g. due 
to the external random perturbations of the system's parameters. The model 
contains only one relaxation rate, however, the power spectral density of the 
current of particles reveals an exact 1/f spectrum in any desirable wide range 
of frequency and can be expressed by the Hooge formula. It is likely that the 
analysis and the generalizations of the model can strongly influence on the 
understanding of the origin of 1 /f noise. 
PACS number(s): 05.40. +j, 02.50.-r, 72.70.+m 
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The puzzle of the origin and omnipresence of 1/f noise - also known as 'flicker' or 'pink' 
noise - is one of the oldest unsolved problem of the contemporary physics. Since the first 
observation of the flicker noise in the current of electron tube more than 70 years ago by 
Johnson fluctuations of signals and physical variables exhibiting behavior characterized 
by a power spectral density S (/) diverging at low frequencies like 1//"^ (5 ~ 1) have been 
discovered in large diversity of uncorrelated systems. We can mention here processes in 
condensed matter, traffic flow, quasar emissions, music, biological, evolution and artiflcial 
systems, human cognition and even distribution of prime numbers (see and references 
herein) . 

1/f noise is an intermediate between the well understood white noise with no correla- 
tion in time and the random walk (Brownian motion) noise with no correlation between 
increments. Widespread occurrence of signals exhibiting power spectral density with 1/f 
behavior suggests that a general mathematical explanation of such an effect might exist. 
However, except some formal mathematical speculations like "fractional Brownian motion" 
or half-integral of a white noise signal [^] no generally recognized mathematical explanation 
of the ubiquity of 1/f noise is still proposed. Physical models of 1/f noise in some physical 
systems are usually very specialized, complicated (see 0-0 and references herein) and they 
do not explain the omnipresence of the processes with 1/ spectrum . 

Note also some mathematical algorithms and models of the generation of the processes 
with 1/f noise ||T0|-[12|. These models also expose some shortcomings: they are very specific, 
formal or unphysical, they can not, as a rule, be solved analytically and they do not reveal 
the origin as well as the necessary and sufficient conditions for the appearance of 1/f type 
fluctuations. 

History of the progress in different areas of physics indicates to the crucial influence 
of simple models on the advancement of the understanding of the main points of the new 
phenomena. We note here only the decisive influence of the Bohr model of hydrogen atom 
on the development of the quantum theory, the role of the Lorenz model as well as the 
logistic and standard (Chirikov) maps for understanding of the deterministic chaos and the 
quantum kicked rotator for the revealing the quantum localization of classical chaos. On 
the contrary, the simple model of 1/f noise, to the best of our knowledge, is absent until 
now. This is one of the reasons of the incomprehensibility of the origin and nature of this 
phenomenon. 

It is the purpose of this letter to present the simplest analytically solvable model of 1/f 
noise which can essentially influence on the understanding of the origin, main properties 
and parameter dependences of the intensity and domination of the flicker noise. Our model 
is a result of the search of necessary and sufficient conditions for the appearance of 1/f 
fluctuations in simple systems affected by the random external perturbations initiated in 
and originated from the observation of the transition from chaotic to nonchaotic behavior 



in the ensemble of randomly driven systems [0]. Contrary to the McWhorter model |T5 



based on the superposition of large number of Lorentzian spectra and requiring a very wide 
distribution of relaxation times, our model contains only one relaxation rate and can have 
an exact 1/f spectrum in any desirable wide range of frequency. 

The simplest version of our model consists of one particle moving in the closed contour. 
The period of the drift of the particle round the contour fluctuates (due to the external ran- 
dom perturbations of the system's parameters) about some average value r. So, a sequence 
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of the transit times tk when the particle crosses some point of the contour is described 
by the recurrent equations 



tk — tfc-l + Tk, 

Tk = Tk-l - 7 (Tfe_i - r) + (7£jfc. 



(1) 



Here 7 <S 1 is the period's relaxation rate, {e^} denotes the sequence of uncorrelated 
normally distributed random variables with zero expectation and unit variance (the white 
noise source) and a is the standard deviation of white noise. 

The intensity of the current of the particle through the section of the contour Lm is 
expressed as 

i{t)^j:s{t-t,) (2) 



where S {t) is the Dirac delta function. The power spectral density of the current (2) is 
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(3) 



where T is the whole observation time interval, fcmin and k 



k q — fc^ii^ k I 

are minimal and maximal 
values of index k in the interval of observation and the brackets (...) denote the averaging 
over realizations of the process. 

Prom Eqs. (1) it follows an expression for the period 



r, = r + (ro-r)(l-7)' + a$:(l-7)'-^'£. 



(4) 



where Tq is the initial period. 

After some algebra we also can easily obtain an explicit expression for the transit times 

tki 



t, = t, + kr+^-Y.[l-{l-ir'-' 
7 1=1 



(5) 



Here is some constant for k ^ 7"^ or tq — r. In the later case is the initial time tg. 
At A; ^ 7~^ Eq. (5) generates the stationary time series and the transit times tk+g and tk 
difference in Eq. (3) is 

( k k+q ^ 

tk^q -h^qr + ^{[l-{l- E (1 - 7)'^'-' + E [1 - (1 - 7)'^"^'-^ ei\,q>0. 

7 [ 1=1 i=k+i } 

(6) 

Substitution of Eq. (6) into Eq. (3) and averaging over realizations of the process or over 
the normal distribution of the random variables £/ yield 



s(f)= lim ±yye'^^f^i-^''''f'<^^i^ 



where 



Y 



[1 - (1 - 7)f E (1 - + E [1 - (1 - 7)'] H , q>0. 



(7) 



(8) 



1=1 
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Summations in Eq. (8) for /c ^ 7 ^ result in 

2 f (1 - 7) [1- (1- 7)1 

7n i-(i-7) 

Expansion of expression (9) in powers of 75 <C 1 is 



s(?) = V- ^ + ^9' + .... <;>o. (10) 

7 O Z 



Note the parity of the function g{q), g{—q) = g{(l), following from Eqs. (6)-(8) at A; — \q\ ^ 

For f^fr={ 
by the integration 



For f <^ fr = (27rr) and f < f2 — I-^I'ko we can replace the summation in Eq. (7) 



S if) =21 j e^2Tr/q-7r2aV^9(9)(^^. (H) 
—00 

Here / = lim (/Cmax — ^min + ^) /T = is the averaged current. Furthermore, at / » 

T— >oo 

/i = 7^/^/7r(7 it is sufficient to take into account only the first term of expansion (10). 
Integration in Eq. (11) hence yields to 1/f spectrum 

S{f) = pY^ /l</</2,/r (12) 

where an is a dimensionless constant (the Hooge parameter) 

2^e-^^ K=^. (13) 
^/^^ a 

The model containing only one relaxation time 7"^ for sufficiently small parameter 7 
can, therefore, produce an exact 1/f-like spectrum in any desirable wide range of frequency. 
Furthermore, due to the contribution to the transit times tk oi the large number of the 
random variables (l = 1,2, ...A;), our model represents a 'long-memory' random process. 
As a result of the nonzero relaxation rate (7 7^ 0) and, consequently, due to the finite 
dispersion of the r period, cx^ = (r|) — (rfc)^ = cr^/27(l —7/2) (2A;7 ^ 1), the model, 
however, is free from the unphysical divergency of the spectrum at / — > 0. So, using an 
expansion of expression (9) at 75 » 1, g{q) = 2q/^'^, we obtain from Eq. (11) the spectrum 
density S{f) = P {2ayT^') for / « /i, /o = r^'/na^ 

Eqs. (11)-(13) describe quite well the power spectrum of the random process (1). As 
an illustrative example in Fig. 1 the numerically calculated power spectral density averaged 
over five realizations of the process (1) is compared with the analytical calculations according 
to Eqs. (9)-(13). Note, that for the used in calculations parameters r, cr, and 7 according to 
Eqs. (11)-(13) we have K = 2, an = 0.04, fi = 5x 10-^ /a = 0.06, fr = 0.08, fo = lx 10-^ 
while spectrum for /i < / < /2 according to Eq. (12) is S{f) = 0.01// and for / << fi, /o 
tends to the constant S{f « fi) = 156. These results are in excellent agreement with the 
numerical analysis. 

This simple, consistent and exactly solvable model can easily be generalized in different 
directions: for large number of particles moving in similar contours with the coherent (iden- 
tical for all particles) or independent (uncorrelated for different particles) fluctuations of the 
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periods, for the non-Gaussian and for the continuous perturbations of the systems' param- 
eters and for some spatially extended systems. So, when an ensemble of particles moves 
in the closed contours and the period of each particle fluctuates independently (due to the 
perturbations by uncorrelated, different for each v particle, sequences of random variables 
{e^}) the power spectral density of the collective current I of all particles can be calculated 
by the above method too and is expressed as the Hooge formula p| 

SU) = I'^- (14) 

Summarizing, simple analytically solvable model of 1/f noise revealing main features and 
parameter dependences of the power spectral density of the noise is proposed. The model 
and its generalizations may essentially influence on the understanding of the origin, the main 
properties and parameter dependences of the intensity and domination of the flicker noise. 

The author acknowledges stimulating discussions with Prof. R. Katilius, Prof. A. Mat- 
ulionis, and Dr. A. Bastys as well as numerical simulation data provided by Mr. T. 
Meskauskas. The research of this publication was made possible in part by support of the 
Alexander von Humboldt Foundation and Lithuanian State Science and Studies Foundation. 
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Captions to the figure of the paper 
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Fig.l. Power spectral density vs frequency of the current generated by Eqs. (1) with 
T = 2, a = 0.04, and 7 = 0.016. Sinuous fine curve is the averaged over five reahzations 
resuh of numerical simulations, the heavy line corresponds to the numerical integration 
of Eq. (11) with g (q) from Eq. (9), and the thin straight line represents the analytical 
spectrum according to Eqs. (12) and (13). 
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Fig. 1 



